Mon. Not. R. Astron. Soc. 000, 1-13 (2010) Printed 28 January 2013 (MN ETgX style file v2.2) 



The reliability of [Cn] as a star formation rate indicator 



Use De Looze 1 , Maarten Baes 1 , George J. Bendo 2 , Luca Cortese 3 and Jacopo Fritz 

1 Sterrenkundig Observatorium, Universiteit Gent, Krijgslaan 281 S9, B-9000 Gent, Belgium 

~UK ALMA Regional Centre Node, Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, 

Un iversity of Manchester, Oxford Road, Manchester M 1 3 9PL, United Kingdom 

- 'European Southern Observatory;Karl-Schwarzschild Str. 2, 85748 Garching bei Muenchen, Germany 



Received 25 December 2010 



ABSTRACT 

The [Cn] 157.74yum line is an important coolant for the neutral interstellar gas. Since [Cn] is 
the brightest spectral line for most galaxies, it is a potentially powerful tracer of star formation 
activity. In this paper we present a calibration of the star formation rate as a function of the 
[Cn] luminosity for a sample of 24 star-forming galaxies in the nearby universe. This sample 
includes objects classified as Hn regions or LINERs, but omits all Seyfert galaxies with a 
significant AGN contribution to the mid-infrared photometry. In order to calibrate the SFR 
against the line luminosity, we rely on both GALEX FUV data, which is an ideal tracer of 
the unobscured star formation, and MIPS 24/rni, to probe the dust-enshrouded fraction of star 
formation. In case of normal star-forming galaxies, the [Cn] luminosity correlates well with 
the star formation rate. However, the extension of this relation to more quiescent (Ha EW 
< 10 A) or ultraluminous galaxies should be handled with caution, since these objects show 
a non-linearity in the LfQij-to-LpiR ratio as a function of Lfir (and thus, their star formation 
activity). 

We provide two possible explanations for the origin of the tight correlation between the 
[Cn] emission and the star formation activity on a global galaxy-scale. A first interpretation 
could be that the [Cn] emission from PDRs arises from the immediate surroundings of star- 
forming regions. Since PDRs are neutral regions of warm dense gas at the boundaries between 
Hn regions and molecular clouds and they provide the bulk of [Cn] emission in most galax- 
ies, we believe that a more or less constant contribution from these outer layers of photon- 
dominated molecular clumps to the [Cn] emission provides a straightforward explanation for 
this close link between the [Cn] luminosity and SFR. Alternatively, we consider the possibil- 
ity that the [Cn] emission is associated to the cold ISM, which advocates an indirect link with 
the star formation activity in a galaxy through the Schmidt law. 

Key words: galaxies: star formation - ISM: lines and bands - infrared: galaxies - ultraviolet: 
galaxies 



1 INTRODUCTION 

The [Cn] line is an important coolant of the neutral interstellar 
medium (ISM), heated through the photoelectric effect on dust 
grains and Polycyclic Aromatic Hydrocarbons (PAHs), that have 
been exposed to ultraviolet (UV) photons. This UV emission can 
originate both in Hn regions, the diffuse neutral and ionized in- 
terstellar medium (Nakagawa et al. 1998; Wolfire 1995) and in 
PhotoDissociation Regions (PDRs) (Crawford et al. 1985; Tielens 
& Hollenbach 1985; Wolfire et al. 1989; Hollenbach et al. 1991; 
Stacey et al. 1991; Bakes & Tielens 1998). Since the [Cn] line is 
generally a very strong line in all star-forming galaxies and in par- 
ticular for low-metallicity objects (Madden 2000), it is a potentially 
powerful star formation rate (SFR) indicator. In this paper we aim 
to quantify the reliability of [Cn] as a star formation rate indicator 
from a sample of galaxies in the nearby universe. 



Despite the strength of the [Cn] line emission, several issues 
might obstruct a direct link with the star-forming activity in galax- 
ies. For instance, the ambiguity concerning the origin of the dom- 
inant heating source for the interstellar gas, contaminates the cor- 
relation between the [Cn] line emission and the SFR. Furthermore, 
this correlation is also affected by the saturation of the [Cn] excita- 
tion at high temperatures and high densities (Kaufman et al. 1999). 
Nevertheless, using [Cn] as a star formation indicator does benefit 
from the fact that it is unaffected by extinction in most cases. Ex- 
ceptions might be extreme starbursts (Luhman et al. 1998; Helou 
2000) or edge-on galaxies (Heiles 1994). Recently, Papadopoulos 
et al. (2010) claimed that the deficiency of [Cn] in Arp 220 is due 
to a high FIR/submm dust optical depth effect and is not caused by 
saturation of the [Cn] excitation in higher density PDRs. Rangwala 
et al. (2010) (private communication) also confirms that [Cn] might 
be affected by higher optical depths in Arp 220, but they found 
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much lower values than Papadopoulos et al. (2010). A thorough 
analysis of this dust obscuration effect at FIR/submm wavelengths 
in other galaxies with a [Cn] deficiency is necessary to give insight 
on this debated issue. 

The Herschel satellite (Pilbratt et al. 2010) is currently mak- 
ing [Cn] observations, with the aim of studying the ISM of nearby 
galaxies. The advent of ALMA will make this line observable at 
high redshifts (for redshift z>2, the [Cn] 157.74 yum line shifts to at- 
mospheric windows at submm/mm wavelengths, which are observ- 
able from ground-based facilities). Stark (1997) and Nagamine et 
al. (2006) already probed the detectability of [Cn] in high-redshift 
galaxies. According to Nagamine et al. (2006) very bright Ly- 
man break galaxies and distant red galaxies are likely candidates 
to have the most prominent [Cn] emission at high redshift. Some 
detections of [Cn] in high-redshift objects already have been re- 
ported (Maiolino et al. 2005; Iono et al. 2006; Maiolino et al. 2009; 
Hailey-Dunsheath et al. 2010; Ivison et al. 2010; Wagg et al. 2010; 
Staceyetal. 2010). 

The [Cn] line has been used previously to probe the physi- 
cal processes in PDRs, but also to detect star formation in nearby 
galaxies (Genzel & Cesarsky 2000; Malhotra et al. 2001; Boselli 
et al. 2002; Luhman et al. 2003; Pierini et al. 2003). Stacey et al. 
(1991) and Pierini et al. (1999) even considered the use of [Cn] as a 
diagnostic for the SFR in non-starburst galaxies. Also Leech et al. 
(1999) found a correlation between [Cn] emission and the star for- 
mation activity in an indirect way. They inferred a direct correlation 
of the [Cn]-to-FIR flux ratio and [Cn]-to-K'-band flux ratio with 
the galaxy lateness, which is correlated, in its turn, with massive 
star formation activity. Boselli et al. (2002) extended this analysis 
and calibrated the [Cn] line flux as a star formation tracer based on 
Hq-+[Nii] line fluxes. They found that the [Cn] line intensity is not 
simply proportional to the star formation rate. Moreover, the scatter 
around the SFR calibration is considerable (an uncertainty of a fac- 
tor ~ 10 when estimating the SFR from the [Cn] line emission). 
Considering the recent improvement in techniques to deal with 
[Nn] contamination (Decarli et al. 2007) and attenuation correc- 
tion (Boselli et al. 2009) of Ha data, we will extend the work done 
by Boselli et al. (2002) and provide a new calibration for the SFR 
relation based on other reliable star formation indicators. The sam- 
ple in Boselli et al. (2002) consisted of 22 nearby late-type galaxies 
spanning a far-IR luminosity range of 10 s < Lfir < 10 105 L Q . We 
will extend this sample to more luminous galaxies up to L FIR ~ 
10 : ' /... 

In case of normal star-forming late-type galaxies, a combina- 
tion of indicators that trace the dust-enshrouded and unobscured 
star formation are able to give a complete picture of the star forma- 
tion activity. Since UV radiation is mainly originating in massive 
young stars, it is an ideal tracer for the unobscured star formation. 
Solely relying on this unobscured fraction will underestimate the 
true star formation activity, as starburst regions are often affected 
by strong attenuation from surrounding interstellar dust clouds. To 
acquire a complete picture of the rate at which stars are formed, 
the dust-enshrouded star formation must be traced as well. The to- 
tal amount of enshrouded star formation can be traced by the to- 
tal IR-luminosity. But also several monochromatic SFR indicators 
have proven to be reliable estimators for the star formation activity. 
In particular, the mid-IR emission at 24/jm shows a general cor- 
relation with the star formation rate (Boselli, Lequeux & Gavazzi 
2004; Calzetti et al. 2005; Wu et al. 2005; Alonso-Herrero et al. 
2006; Perez-Gonzalez et al. 2006; Relano et al. 2007; Calzetti et 
al. 2007; Bavouzet et al. 2008; Zhu et al. 2008; Rieke et al. 2009; 
Salim et al. 2009; Kennicutt et al. 2009). 



In this paper, we investigate the correlation between [Cn] lu- 
minosities and the SFR, estimated from the GALEX FUV and 
MIPS 24 ytim or a combination of these SFR tracers. Based on the 
strength of these correlations, we will choose the SFR indicator 
showing the tightest correlation with [Cn], which will serve as a 
reference SFR diagnostic in our analysis to calibrate the SFR-[Cn] 
relation. §2 describes the sample selection and data acquisition. In 
§3, we estimate the SFR either from the monochromatic 24yum lu- 
minosity or in combination with the GALEX FUV luminosity. Fi- 
nally, §4 discusses the reliability and applicability of [Cn] as a star 
formation rate indicator and also addresses the nature of [Cn] emis- 
sion in galaxies. §5 briefly summarizes our main conclusions. 



2 SAMPLE AND DATA 

2.1 Sample selection 

The selection is based on the galaxy sample in Brauher et al. (2008), 
who assembled all galaxies with available [Cn] data from the ISO 
archive. More specifically, we restrict our sample to those galax- 
ies for which the [Cn] line fluxes do not correspond to upper limits 
(i.e. non-detections) and which are classified in this paper as unre- 
solved in the far-IR with respect to the ~ 75" ISO LWS beam. This 
latter criterion implies that an aperture correction is not required for 
the [Cn] line flux. Although a substantial fraction of the flux from 
an on-axis point source is diffracted out of the aperture beyond the 
diffraction limit at about 1 10 fim for all LWS detectors, these losses 
will be cancelled out in the calibration process, provided it is ap- 
plied to point sources observed on-axis (Gry et al. 2003). 

From this sample of 153 unresolved galaxies, we retain all 
galaxies that have been observed in both GALEX FUV and MIPS 
24yum bands, which gives us a final sample of 39 unresolved galax- 
ies. 

Some galaxies have multiple [Cn] data, corresponding to aper- 
tures taken at different positions. For every galaxy, we only include 
one data point, preventing the introduction of any bias in our sam- 
ple. If the different apertures show some overlap in area, we choose 
the aperture corresponding to the [Cn] flux with the smallest uncer- 
tainty. If the overlap is insignificantly small, we add all fluxes from 
the apertures taken at different positions. The fluxes from the differ- 
ent apertures for the GALEX FUV and MIPS 24 jim are summed 
accordingly. This kind of flux summation has only been applied for 
two galaxies in our sample: NGC4651 and NGC7217. 

Since FUV data are affected by galactic extinction, we have 
to apply appropriate corrections. At a galactic latitude b = 11. °2 
(Clements 1983), NGC 1569 suffers from a large amount of Galac- 
tic extinction. Inconsistencies among the different reported esti- 
mates for the V-band extinction (e.g. Ay = 1.79 (Israel 1988); A\ = 
1.61 (Devost et al. 1997); A v = 2.32 (Schlegel et al. 1998)), made 
us decide to remove NGC 1569 from our sample, considering we 
want to avoid the introduction of any bias. This reduces our sample 
to 38 galaxies. 

2.2 GALEX data 

Far-ultraviolet (FUV, /1=1539A, A/1=442A) observations of galax- 
ies in our sample have been obtained from the GALEX GR4/5 pub- 
lic release. In order to derive accurate FUV photometry we only 
used fields with integration times greater then 800 sec, obtained 
as part of the Nearby Galaxy Survey, the Medium Imaging Survey 
or Guest Investigator programs (70 galaxies from our unresolved 
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sample of 153 galaxies have been observed by GALEX). We per- 
formed aperture photometry on the intensity maps produced by the 
standard GALEX pipeline. Details about the GALEX pipeline can 
be found in Morrissey et al. (2007). The GALEX FUV data have 
been corrected for galactic extinction according to Schlegel et al. 
(1998), using the extinction relation obtained from Cardelli et al. 
(1989). 

2.3 MIPS 24pm data 

The 24 micron images were created from raw data produced by 
the Multiband Imaging Photometer for Spitzer (MIPS; Rieke et al. 
2004) on the Spitzer Space Telescope (Werner et al. 2004). The raw 
data were taken by a variety of programs using either the photome- 
try mode (which produces 5 times 5 arcmin maps) or the scan map 
mode (which produces images that are typically 1 degree in length). 
We used the MIPS Data Analysis Tools (Gordon et al. 2005) along 
with additional processing steps to produce the final images. The 
individual data frames were processed through droop correction 
(to remove excess signal in each pixel proportional to the signal in 
the entire array), non-linearity correction, dark current subtraction, 
scan mirror position-dependent and position-independent flatfield- 
ing, and latent image removal steps. Background emission from 
the zodiacal light and additional scattered light that is related to the 
scan mirror position is then subtracted out. For objects observed 
in the scan map mode, the data from separate astronomical ob- 
servation requests were used to make mosaics taken at individual 
epochs. These images were then subtracted from each other to iden- 
tify asteroids, and the regions with asteroids were then masked out. 
Preliminary mosaics of each object were made using all data for 
each object to identify pixels in individual frames with values that 
are statistical outliers compared to other co-spatial pixels in other 
frames. These pixels were masked out, and then final mosaics were 
made using an image scale of 1 .5 arcsec pixeL 1 . Any residual back- 
ground in the images was measured in multiple regions surround- 
ing the optical disc and was then subtracted from the image. The 
final images are then calibrated using the conversion factor given 
by Engelbracht et al. (2007), which has an uncertainty of 4%. The 
FWHM of the PSF is 6 arcsec (Spitzer Science Center 2007). 

2.4 Flux determination 

We determine GALEX FUV and MIPS 24 pm fluxes within the 
same aperture as the ~ 75" ISO LWS beam. Variation in the 
shape of the point spread function used by the different instruments 
should not be an issue, as the necessary corrections have been ap- 
plied through the calibration on a point source. 

The GALEX FUV and MIPS 24pm fluxes were determined 
using the source-extracting code SExtractor (Berlin & Arnouts 
1996). To align the apertures with the exact same position at which 
the [Cn] fluxes were determined, we have created a fake object in 
IRAF (using the task mkobjects) at the same position coordinates 
of the actual pointing of the [Cn] observation. Feeding SExtrac- 
tor with this artificially created map and the FUV or 24pm image, 
assures us that the FUV and 24pm fluxes are extracted within the 
same area on the sky as the [Cn] fluxes. 

2.5 Spectral classification 

Considering we will calibrate the [Cn] luminosity against the SFR, 
derived either from the monochromatic 24pm flux or in combi- 



□ Hll/storburst 

LINER 
a Seyfert 




L°9[L[ai] (erg S'')} 



Figure 1. L[Cu] as a function of L^um^ Upper panel: the three subsamples 
are represented (blue squares = Hn/starburst, cyan diamond = LINER, red 
triangle = Syl/Sy2). For reasons of clarity, we have omitted the errorbars 
from this plot. The best fitting line through the data points are indicated with 
a blue dotted, cyan dashed and red dashed-dotted line for the Hn/starburst, 
LINER and Seyfert subsample, respectively. The mean trend for the com- 
plete sample is represented by a black solid line. The data points have a 
dispersion of 0.5 1 dex around this mean trend. Lower panel: the reduced 
sample, combining Hn, starburst and LINER galaxies. The dotted line rep- 
resents the best fitting line through the data points with a dispersion of 0.3 1 
dex. 



nation with the FUV flux, we have to make sure that the star- 
forming regions are the dominant contributor to this MIR emis- 
sion. Therefore, we will classify our sample according to their 
nuclear spectral classification (see Table 1). We have made three 
different subdivisions. The first subgroup encompasses all objects 
resembling features typical of Hu regions or starbursts. The sec- 
ond and third subset will include LINER and Seyfert galaxies, 
respectively. For this spectroscopic classification, we mainly rely 
on Veron-Cetty & Veron (2010), who have composed a compila- 
tion of AGN host galaxies. For the identification of the remaining 
sources and the objects with an uncertain classification in Veron- 
Cetty & Veron (2010), we have adopted the classification based 
on the optical spectra in Veilleux et al. (1995) and Ho, Filippenko 
& Sargent (1997), who both use similar selection criteria in their 
classification procedure. Additionally, for some sources we used 
the classification from Veron-Cetty & Veron (1986). Considering 
that the optical spectroscopic classification for FIR-luminous ob- 
jects not always offers a clear distinction between starburst- or 
AGN-dominated nuclear activity, we additionally examine spec- 
troscopic L-band data (Risaliti 2004) and radio properties (Baan 
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& Klockner 2006), if available for those objects. Based on sev- 
eral L-band diagnostics, Risaliti (2004) found a dominant AGN- 
contribution to the energy output in the infrared for IRAS 19254- 
7245, IRAS 2055 1-4250 and one of the two optical nuclei in IRAS 
23128-5919. Baan & Klockner (2006) found from their radio ob- 
servations that an AGN is the dominant power source in the nucleus 
of CGCG1510. 8+0725. Four sources remain unclassified, but since 
they all host Hn regions (Cartwheel: Gao et al. 2003, NGC1317: 
Crocker et al. 1996, NGC4189 and NGC4299: Hodge & Kenni- 
cutt 1983), we assign them to the Hn/starburst group. 

Our sample of 38 galaxies contains 16 galaxies for which 
Hn regions or starbursts dominate their central regions, while the 
remaining sources in our sample are host galaxies of an AGN 
(14 galaxies) or LINER (8 galaxies). For these different nuclear 
regimes, we will examine the power sources that mainly contribute 
to the 24 jim emission. Plotting L[ Cll ] as a function of the MIPS 
24/im luminosity for the three different subsets, we immediately 
deduce a significant contribution from the AGN in Seyfert galaxies 
to the 24pm emission (see the upper panels in Figure 1). The dis- 
persion around the mean trend for the complete sample of galaxies 
is 0.51 dex in Figure 1. For galaxies classified as a LINER, star- 
bursts appear to be the main power supply for the MIPS 24pm. 
Although the spectral line emission in some LINERs might also 
be generated by a quiescent AGN, this contribution does not ap- 
pear to be significant for the LINER galaxies in our sample. Be- 
cause of the insignificant AGN contribution to the MIR emission in 
LINERs, we perform our calibration analysis on a reduced sample 
that combines the two subsets (Hn regions, starbursts and LINERs), 
which all have a dominant contribution from starburst to the MIR 
emission. This final sample consists of 24 galaxies of which 16 
show features of Hn regions and 8 other galaxies are classified as 
a LINER. Plotting again the same correlations for this final sam- 
ple (see Figure 1, bottom panel), we can quantify that most of the 
scatter in previous plots was due to the contribution of an AGN to 
the 24pm emission. Indeed, the dispersion around the mean trend 
reduces to 0.31 dex in Figure 1 (bottom panel). 

2.6 Sample description 

The galaxies in our final sample span a range of almost 4 orders of 
magnitude in total infrared luminosity (L T1R ), from ~ 1.2 x 10 42 to 
~ 2.7 x 10 45 erg s _1 (Ltir was calculated based on equation (5) in 
Dale & Helou 2002, relying on the IRAS 25, 60 and 100 pm fluxes). 
In this final sample of 24 galaxies (Hn regions, starbursts and LIN- 
ERs), five galaxies are classified as LIRG (Ltir > 10" L Q ). Figure 
2 demonstrates this wide range in Ltir and shows its dependence 
on the / v (60pm)-to-/ y (100pm) ratio, which is a proxy for the ef- 
fective dust temperature. As most galaxies in our sample have high 
effective temperatures, the IRAS 100pm flux samples the SED peak 
of dust emission rather well. This was also confirmed by Kenni- 
cutt et al. (2009), as they found that the IRAS fluxes give a reliable 
representation of the total infrared luminosity for galaxies obeying 
/ y (60pm)// v (100pm) > 0.4 . Only six galaxies in our sample have a 
/,,(60pm)-to-/ v (100pm) ratio below this value. 

Our sample covers a large range in distance, going from 3.9 to 
139.1 Mpc. Four galaxies (NGC4189, NGC4293, NGC4299 and 
NGC 7714) are in common with the sample that was used in Boselli 
et al. (2002). 

Table 1 gives an overview of all relevant information for our 
sample galaxies. We arranged their properties as follow: 
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Figure 2. The ratio / v (60//m)-to-/,,(100/im), which is a proxy for the effec- 
tive dust temperature, as a function of the total infrared luminosity Ltir. 

• Column 1 : galaxy name, NGC, IC, CGCG, UGC or IRAS, for 
the sample galaxies. 

• Column 2: distance [Mpc], derived from the Nearby Galaxies 
Catalogue (Tully 1988) and for other galaxies, from their recession 
velocity (from NED) assuming H = 70 km s -1 Mpc -1 . 

• Column 3: morphological type, adopted from de Vaucouleurs 
etal. (1991). 

• Column 4, 5: classification of the nuclear spectrum: Hn or star- 
burst, LINER, Seyfert (1 or 2), and the corresponding references for 
this classification. 

• Column 6: / v (60/im)-to-/ v (100pm) ratio of the IRAS fluxes, a 
proxy for the effective dust temperature. 

• Column 7: Total-IR luminosity [erg/s] for the whole galaxy, as 
calculated from equation (5) in Dale & Helou (2002). 

Table 2 summarizes all relevant published data and measured 
fluxes within the ~ 75" ISO LWS beam for the analysis in this 
paper. Summing up: 

• Column 1 : galaxy name, NGC, IC, CGCG, UGC or IRAS, for 
the sample galaxies. 

• Column 2, 3: Actual pointing position: right ascension (RA) 
and declination (DEC) from the [Cn] observation, both in decimal 
degrees (from Brauher et al. 2008). 

• Column 4: [Cn] line flux, in [fW/m 2 ], and the uncertainty on 
this flux measurement (from Brauher et al. 2008). 

• Column 5: GALEX FUV flux, [mjy], and the corresponding 
uncertainty. Both values have been corrected for galactic extinction, 
according to Schlegel et al. (1998) and Cardelli et al. (1989). 

• Column 6: Spitzer MIPS 24 pm flux, [mJy], and the corre- 
sponding uncertainty, including both uncertainties in the flux ex- 
traction (< 2%) and the calibration (4%, Engelbracht et al. 2007). 

Although this final sample of 24 galaxies, gathering objects 
classified as Hn region, starburst or LINER, is neither statistically 
significant nor representative for the whole of galaxies with diverg- 
ing properties, this sample enables us to make a first preliminary 
analysis of the diagnostic capabilities of [Cn] in tracing the star 
formation activity in galaxies. For a thorough analysis based on a 
more extensive sample of galaxies, we will have to wait until the 
completion of all [Cn] surveys with the Herschel Space Observa- 
tory. Nevertheless, the broad range of optical and infrared lumi- 
nosities covered in our sample enables us to infer the reliability of 
[Cn] as a star formation indicator in star-forming galaxies spanning 
almost 4 orders of magnitude. 
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Table 1. Properties for the galaxies in our sample. 



Name Distance Type 1 Spectral type Ref 2 IRAS 60 / IRAS 100 logL T lR 

(Mpc) log (erg/s) 

(1) (2) (3) (4) (5) (6) (7) 



Cartwheel 


129.3 


RING 


Hn/sb 


3 

no 


0.46 


44.30 


NGC 0520 


27.8 


Pec 


Hn/sb 


3 


0.66 


44.41 


NGC 0625 


3.9 


SB(s)m? edge-on 


Hn/sb 


1 


0.56 


42.08 


NGC 0660 


11.8 


SB(s)a pec 


LINER 


1, 2, 3 


0.64 


44.03 


NGC 0695 


139.1 


SO? pec 


Hn/sb 


2 


0.58 


45.28 


NGC 0986 


23.2 


SB(rs ab 


Hn/sb 


4 


0.49 


44.32 


UGC 02238 


92.3 


Im? 


LINER 


2 


0.54 


44.92 


NGC 1156 


6.4 


IB(s)m 


Hn/sb 


3 


0.50 


42.48 


NGC 1266 


31.3 


(R') SB(rs) 0a0 pec: 


LINER 


1.2 


0.81 


44.07 


NGC 1275 


75.2 


Pec 


Seyfert 


1, 3 


1.02 


44.49 


NGC 1317 


16.9 


SAB(r)a 


Hn/sb 


no 


0.34 


43.28 


NGC 1569 


1.6 


IBm 


Hn/sb 


3 


0.98 


42.15 


IRAS 05189-2524 


182.3 




Seyfert 


1. 2 


1.20 


45.74 


UGC 03426 


57.9 


SO: 


Seyfert 


1 


1.12 


44.58 


NGC 2388 


59.1 


S? 


Hn/sb 


2 


0.67 


44.83 


NGC 4041 


22.7 


SA(rs)bc: 


Hn/sb 


3 


0.41 


44.06 


NGC 41 89 


16.8 


SAB(rs) cd? 


Hn/sb 


no 


0.34 


43.24 


NGC 4278 


9.7 


El-2 


LINER 


3 


0.36 


42.09 


NGC 4293 


17.0 


(R) SB(s) 0/a 


LINER 


3 


0.44 


43.32 


NGC 4299 


16.8 


SAB(s)dm: 


Hn/sb 


no 


0.33 


43.17 


NGC 4490 


7.8 


SB(s)d pec 


Hn/sb 


3 


0.60 


43.51 


NGC 4651 


16.8 


SA(rs)c 


LINER 


3 


0.35 


43.44 


NGC 4698 


16.8 


SA(s)ab 


Seyfert 


1,3 


0.31 


42.65 


IC4329A 


68.8 


SA0A+: edge-on 


Seyfert 


1 


1.22 


44.59 


NGC 57 13 


30.4 


SAB(rs)bc pec 


Hn/sb 


4 


0.57 


44.43 


CGCG 1510.8+0725 


55.7 




Seyfert 


6 


0.66 


44.76 


NGC 6221 


19.4 


SB(s)c 


Seyfert 


1 


0.50 


44.37 


NGC 6240 


104.8 


10: pec 


LINER 


1. 2 


0.82 


45.43 


IRAS 19254-7245 


264.3 




Seyfert 


1,5 


0.95 


45.68 


NGC 6810 


25.3 


SA(s)ab: 


Seyfert 


4 


0.52 


44.28 


IRAS 20551-4250 


184.1 




Seyfert 


1,5 4 


1.29 


45.55 


NGC 7217 


16.0 


(R)SA(r)ab 


LINER 


1,3 


0.27 


43.46 


NGC 7469 


69.9 


(R')SAB(rs)a 


Seyfert 


1, 2 


0.74 


45.25 


IRAS 23128-5919 


191.0 




Seyfert 


1.5 s 


0.98 


45.58 


NGC 7552 


19.5 


(R')SB(s)ab 


Hn/sb 


1 


0.71 


44.54 


NGC 7582 


17.6 


(R')SB(s)ab 


Seyfert 


1 


0.67 


44.25 


NGC 7714 


36.9 


SB(s)b: pec 


Hn/sb 


1. 2 


0.90 


44.31 


IRASF23365+3604 


276.2 




Seyfert 


1 


0.84 


45.74 


NGC 7771 


61.1 


SB(s)b: pec 


Hn/sb 


2 


0.49 


45.00 



1 The galaxy type has been adopted from de Vaucouleurs et al. (1991). If no galaxy type is mentioned, there was no classification available in the literature. 

2 The references of the spectral type: 1: Veron-Cetty & Veron (2010); 2: Veilleux et al. (1995); 3: Ho, Filippenko & Sargent (1997); 4: Veron-Cetty & 
Veron (1986); 5: Risaliti (2004) 6: Baan & Klockner (2006) 

3 "no" indicates that no classification was found in the literature. These objects were classified in the Hn/starburst group based on detected Hn regions in 
their nuclei. Section 2.5 discusses the classification for these galaxies more into detail. 

4 Based on optical spectra IRAS 2055 1-4250 was classified as a starburst galaxy (Veron-Cetty & Veron 2010), but according to L-band diagnostics a 
significant AGN contribution is present (Risaliti 2004). 

5 This galaxy has a pair of optical nuclei, which is probably a remnant of a merging process. Veron-Cetty & Veron (2010) identified this galaxy as Hn 
region, while Risaliti (2004) reported that one nucleus shows AGN features and the other nucleus resembles a Hn region. Since both nuclei fall within the 
ISO LWS beam, we have classified this object as a Seyfert galaxy, based on the AGN-contribution from at least one nucleus. 



3 RELIABLE REFERENCE SFR TRACER 

In order to compare the different star formation relations, it is im- 
portant that they are calibrated in the same way. Therefore we will 
relate all relations to the same initial mass function (IMF). This 
reference IMF will be the Kroupa (2001) IMF characterized by a 
power law £(m) oc m " with slope a = 2.3 for 0.5-100 M Q and a = 



1.3 for 0.1-0.5 M Q . We prefer the Kroupa (2001) IMF rather than 
the more commonly used Salpeter (1955) IMF (with a single slope 
a = 2.35 for 0.1-100 M Q ), since the former better correlates with 
IMF observations in the Galactic field (Chabrier 2003; Kroupa & 
Weidner 2003). 

Analyzing the correlations between L[c„] and the SFR requires 
a complete threefold approach: 
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Table 2. Sample galaxies: data 



Name RA DEC [Cn] GALEXFUV MIPS 24 

(deg) (deg) (fW/m 2 ) (mjy) (mjy) 

(1) (2) (3) (4) (5) (6) 



\ d.1 LW11CC1 


9 4171 


-33 7239 


o 


j 5 




o 


.02 


i q«n + n 009 


63 7 + 2 7 


Nnr 0590 


21 1458 


+3 7950 




83 







.07 


7S7 + 009 

rJ.IOI ZL U.UU7 


71 1Q 1 + 7 
J7.J zn oj./ 


Nnr 0695 


23 7799 


-41 4357 


j 


.00 




o 


.08 


o + n 01 fi 
o.Ajy ZL U.U1U 


776 0+117 

/ / U.U ZZZ 1 ._ 


Nnr 0660 


25 7595 


+ 1 1 64S8 

TlJ .UlJO 


7 


.67 




o 


.19 


186 + 007 


9Q«Q Q + 1 1 S 7 

i007.7 ZL 1 1 J. / 


Nnr 0695 


27.8091 


+22 5825 


9 


.04 




o 


.08 


3 oao + n 003 


672 4 + 27 


Nnr oo86 

L \ V J V. U7UU 


38 3922 


-39 0462 


3 


.04 




o 




064 4 01 1 


9155 + 04 4 

ZJJJ.7 X 7t,t 


unr 02238 


41.5726 


+ 13 0958 

l 1 J.W7JO 


2 


.35 




o 


.09 


448 + 017 


406 7 + 200 


Nnr 1 1 56 


44 927 1 


+75 7175 

T^Z. J .Z. J / J 


j 


.56 




o 


.14 


ni. ^07 4 1 1 S 


447 1 4 1 8 

'11 1.1 i 1 OA/ 


Nnr 1 966 

\W}K , 1ZOO 


4Q 0019 


-2 4271 


o 


so 




I) 


1 


916 + 00Q 

U.ZIO X U.UU7 


804 1 + 19 1 


Nnr i ii 5 


49 9525 


+41 5105 

Tt 1 .J 1 UJ 


y 


.16 




o 


.07 


4 7 Aft + 014 


9814 7 + 119 7 

Z.OlT^./ i II-./ 


Nnr 1117 


SO 6869 


-37 1027 


o 


Q 1 







05 


1 .7 / 4 X U.ul J 


907 9 + 8 4 
ZU / .Z. X 0.4 


Nnr i 560 


67 7046 


+64 8478 

TUi.Ot / O 


S 


.99 




o 


.16 


105 917 + 114 


A405 8 + 960 


TR AS 051 80-9594 


SO 9558 

OVJ . _ 1 JO 


-25.3624 


o 


14 




o 


0° 


19640 004 

U. 1ZU i V/ . out 


9786 6+1116 

Z. / OvJ. U 31 11 l.V) 


TTnr 01496 


93 9016 


4-71 0176 


o 


41 







.02 


669 4 090 


9918 1 + 80 7 

i.ZJO.1 i 07. / 


Nnr 9188 


1 12 2229 


+ T.T. 81 89 

T J J.O 1 OZ. 


y 


.74 




o 


.08 


1 1S 4 OOS 


1 9 + fti 5 


Nnr 404 1 

1 \ V J v. . tut 1 


1 SO 5501 


469 1 161 

l LIZ. . 1 JUJ 


3 


.48 




o 


.05 


1 880 4 001 


ii 50 7 + 46 5 


Nnr 41 80 


1 81 4465 

1 O J I ' IUJ 


+ 13 4240 


o 


94 




o 


.06 


7 94^ + 00 1 


295 7 4 12 


Nnr 4978 


1 85 0976 

1 OJ .V '_ / u 


49Q 9894 

TZ.^.Z.OZ.'T 


o 


25 




o 


.03 


881 4 010 


49.1 ±21 


Nnr 49 oi 


1 85 1050 


+ 18 3844 


o 


28 




o 


.02 


111 4 001 

U. ill i u.UUJ 


<S90 1 + 90 


Nnr 49QQ 


1 85 41 80 

1 O J .*+ 107 


+ 11 5012 


o 


93 




o 


.09 


QQft -L. 010 
O.77U X U.UJu 


210 4 4 8 5 


Nnr 44Q0 

1> \JV_« 44:7 \J 


1 87 6516 

Jo / .LUjO 


441 6108 

+41 .OJ70 


4 


-3? 







1 1 


1 948 + 099 


1 104 8 + 55 


Nnr 465 i 


190 9271 


+ 16.3942 


2 


.03 




o 


.08 


4 048 4 015 


181 1 + 15 5 

JOJ.l i 1J.J 


Nnr 465 i 


190 9442 


416.3972 


o 


50 







.06 


1 «Q 4. 010 

l.JJJ i W . 1 J 1 w 


79 £ + 3 3 


Nnr 465 i 


190 91 


416 3914 


o 


.76 







.09 


i 719 + o 010 


88 + 1 6 


Nnr 465 i 


1~r*ri 1 frnm 1 i T"\r*ci 1~i on c 

LULal 11U111 till LHJ31L1U113 






.29 







1 } 


7 1 1 4 4 09 1 


550 7 + 99 9 




i vz.uyov 


4o.4o / D 


u 






u 


.\JZ 


u.zoy ± U.UU4 


1/t Q -i- 1 5 


IC4329A 


207.3304 


-30.3095 





.15 


± 





.02 


0.073 ± 0.003 


1945.3 ± 77.9 


NGC5713 


220.0478 


-0.2905 


4 


.68 


+ 





.13 


4.061 ± 0.020 


2144.6 ± 85.9 


CGCG 1510.8+0725 


228.3053 


47.2265 





.33 


+ 





.05 


0.086 ± 0.004 


514.3 ± 20.7 


NGC6221 


253.1913 


-59.2167 


6 


.64 


+ 





.25 


4.527 ± 0.063 


4528.1 ± 181.3 


NGC6240 


253.2450 


+2.4013 


2 


.72 


± 





.06 


0.637 ± 0.004 


2735.4 ± 109.5 


IRAS 19254-7245 


292.8400 


-72.6319 





.26 


+ 





.07 


0.162 ±0.005 


1059.0 ± 42.5 


NGC6810 


295.8915 


-58.6556 


3 


.83 


+ 





.13 


0.553 ± 0.010 


2882.1 ± 115.4 


IRAS 20551-4250 


314.6116 


-42.6518 





.41 


+ 





.03 


0.600 ± 0.006 


1450.9 ± 58.2 


NGC7217 


331.9775 


+31.35833 





.66 


± 





.06 


1.557 ± 0.008 


197.9 ± 8.0 


NGC7217 


331.95833 


+31.35972 





.62 


+ 





.07 


1.331 ± 0.008 


185.2 ± 7.5 


NGC7217 


total from all positions 




1 


.28 


+ 





.09 


2.888 ± 0.012 


383.1 ± 15.5 


NGC7469 


345.81417 


+8.87361 


2 


.27 


± 





.03 


6.963 ± 0.020 


4421.9 ± 177.0 


IRAS 23128-5919 


348.9454 


-59.0544 





.58 


± 





.05 


0.557 ± 0.005 


1242.2 ± 49.8 


NGC7552 


349.045 


-42.5844 


6 


.37 


+ 





.15 


3.746 ±0.011 


9443.8 ± 377.9 


NGC7582 


349.5988 


-42.3703 


4 


.12 


+ 





.13 


0.900 ± 0.006 


5516.4 ± 220.8 


NGC7714 


354.0612 


+2.1550 


1 


.83 


+ 





.10 


7.893 ± 0.018 


2131.7 ± 85.4 


IRASF23365+3604 


354.7554 


+36.3528 





.16 


± 





.02 


0.197 ± 0.006 


585.2 ±23.5 


NGC7771 


357.8534 


+20.1119 


2 


.98 


+ 





.09 


0.906 ±0.012 


1216.5 ±48.8 



• We first analyze the simple relation of the SFR as a function 
of the [Cn] luminosity. Although this approach is most straight- 
forward, the data might be biased towards larger distances. Since 
the distances to the galaxies in our sample range from 3.9 to 139.1 
Mpc, we have to eliminate the distance bias in this luminosity ver- 
sus luminosity relation. 

• Plotting SFR/47r£) 2 as a function of the [Cn] flux, the distance 
bias is removed. If a true [Cn]-SFR correlation is present, it will be 
confirmed in this plot. 

• The dispersion around the mean trend is best quantified in a 
plot of the [Cn] luminosity versus the residual of the best fitting 
relation between L[qh and the SFR. The SFR relation showing the 



tightest correlation with the [Cn] luminosity will be characterized 
by the lowest dispersion in the residual plots. 

3.1 SFR calibrated against TIR 

Taking into account that the 24yum emission from a galaxy better 
correlates with the SFR as opposed to monochromatic FIR mea- 
surements at longer wavelengths (70 and 160yum) (Calzetti et al. 
2010), we have only estimated the SFR, which is used for the fi- 
nal calibration of the SFR-[Cn] relation, from 24yum data (or FUV 
luminosities which have been corrected for attenuation based on 
24yum data) in this work. To obtain a second independent SFR di- 
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agnostic, we have investigated whether TIR luminosities for these 
galaxies could be used as a calibrator. Since we only have IRAS 
fluxes covering the whole galaxy for all sample objects to estimate 
the TIR luminosity, we checked whether an aperture correction was 
required. Although the galaxies in our sample were selected based 
on the criterion that they are unresolved in the FIR with respect to 
the ISO LWS beam, many of them are resolved within the ISO beam 
in optical data. Therefore, we visually inspected the MIPS 24yum 
images and could indeed quantify that also in the MIR a significant 
fraction of the total flux for several galaxies fell outside the ISO 
beam aperture. Subsequently, we verified the availability of MIPS 
70 and 160//m data in the Spitzer archive, to determine the correct 
TIR luminosities which correspond to the portion of the flux that 
falls within the same ISO beam. Considering that only 28 of the 38 
sample galaxies have full coverage in all Spitzer bands, and some 
of those galaxies were only observed in small fields which does not 
allow proper convolution with the appropriate kernel to match the 
PSF's in all Spitzer bands, we decided not to use these Spitzer data. 
Moreover, the varying sizes of the sample galaxies would imply 
a different aperture correction (some are point sources and others 
are clearly more extended), introducing an additional bias to the 
TIR luminosity derived from these Spitzer data. While the MIPS 
24/im clearly demonstrate that the majority of our sample galax- 
ies requires some kind of aperture correction, we have verified that 
the data currently at hand are not able to properly correct for these 
aperture issues. Any empirical correction can significantly bias the 
correlations between the [Cn] luminosity and the SFR derived from 
this aperture-corrected TIR luminosity and does not allow us to 
draw any firm conclusion. Therefore, this paper only investigates 
the correlation between [Cn] and the SFR, when estimated from 
the 24/im or attenuation-corrected FUV flux. 

3.2 SFR calibrated against MIPS 24 yum 

In this section, we predict the star formation rate solely relying on 
the MIR emission at 24/zm, to trace the dust-enshrouded fraction of 
star formation. 

First, we have to choose an appropriate star formation rela- 
tion. Several SFR relations based on 24 /an emission have been 
published. Some authors claim a linear relation between the SFR 
and L 2 4 pm (Wu et al. 2005; Zhu et al. 2008), while others argue 
that only a non-linear relation can account for the increasing 24yum 
emission in more IR-luminous galaxies which form stars at a higher 
rate (Wu et al. 2005; Perez-Gonzalez et al. 2006; Alonso-Herrero 
et al. 2006; Calzetti et al. 2007; Relano et al. 2007; Zhu et al. 2008; 
Rieke et al. 2009; Boquien et al. 2010). As some galaxies in our 
sample are situated at the high luminosity end, we should be aware 
that this non-linear effect might affect our results. We choose to 
estimate the SFR from the relation in Rieke et al. (2009). 

Rieke et al. (2009) use a combination of a linear trend at lower 
TIR-luminosities and a non-linear relation above a threshold of Ltir 
> 10 11 L e or L 2 4 pm > 1.3 x 10 10 L Q . Both relations were derived 
from a sample of dusty, luminous star-forming galaxies and hold 
up to Ltir = 2 x 10 12 L Q , which is higher than the most luminous 
galaxy in our sample (NGC 6240; 7.0 x 10" L Q ). The overlap be- 
tween our range in TIR luminosities and the similar characteristics 
of both samples make us confident that this relation gives a good in- 
dication of the true star formation rate in our sample galaxies. The 
two relations from Rieke et al. (2009) are (assuming our reference 
Kroupa 2001 IMF): 



1 1 1 1 
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Figure 3. Upper panel: SFR 2 v„ (SFR relation from Rieke et al. (2009)) 
as a function of L[Cnj, central panel: SFR24^ m /47rD 2 as a function of [Cn] 
flux, bottom panel: residuals of the SFR^^m-^Qi] plot. In this residual plot, 
Hii/starburst and LINER galaxies are indicated as blue filled squares and red 
diamonds, respectively. 



for 4 x 10 42 < L 24//m < 5 x 10 43 erg s" 1 and 

SFR = 2.04 x 10~ 43 x L 24 ,„„(2.03 x 10- 44 L 24 „„,) 048 (2) 

for L 24 ,, m > 5 x 10 43 erg s -1 , where SFR and L 24 ,, m have units M Q 
yr~' and erg s -1 , respectively. 

After estimating the SFR from these relations, we verify the 
strength of the correlation between L[c„j and the SFR through a 
three-pronged approach. First, we plot the derived SFR as a func- 
tion of L[Qi] (see Figure 3, upper panel). The best fitting line 



SFR = 2.04 x 10~ 43 L 



2Afim 



(i) 
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through the data points corresponds to: 

log SFR 2> „ = 1 .040 log L [C n] - 42.535 (3) 

where the units of SFR^,,, and L[ Cl i] are in M Q yr~' and erg s -1 , 
respectively. The lcr dispersion of the individual galaxies around 
the best fitting line is 0.30 dex. This lcr dispersion has been cal- 
culated as the standard deviation of the logarithmic distance from 
the data points to the best fitting line weighted by the error on the 
individual quantities. Since the observed trend in this plot of lumi- 
nosity versus SFR (and thus indirectly luminosity) might be biased 
by the range in distances among our sample galaxies, we remove 
this distance bias in the central panel of Figure 3. In this figure, we 
plot the SFR-to-47rD 2 ratio as a function of the [Cn] flux. Since this 
distance-independent plot still shows a similar correlation between 
[Cn] and the SFR and the dispersion 0.30 dex is of the same order 
as in the first panel, we safely conclude that the observed correla- 
tion in the upper panel of Figure 3 is not due to a distance bias in 
the data set. The dispersion around the mean trend (see Equation 3) 
is better represented in the bottom panel of Figure 3, which shows 
the residual plot. 

3.3 SFR calibrated against GALEX FUV + MIPS 24 yum 

Before the GALEX FUV flux can be used as a direct quantifier of 
the star formation activity in a galaxy, a correction for internal dust 
attenuation has to be applied. The specific amount of extinction that 
is affecting the UV light is difficult to quantify, because different 
galaxies will be characterized by different extinction laws and a 
variation of stellar and dust distributions. In this section, we will 
rely on SFR relations that use a combination of the GALEX FUV 
flux and the monochromatic MIR emission at 24/jm to obtain a 
complete picture of the star formation activity. 

Two SFR relations have been reported, making use of both 
photometric data. First, Zhu et al. (2008) derived a linear relation to 
estimate the attenuation-corrected FUV-flux based on the observed 
FUV and data: 



iFUV.con U<3] - ^FUVobst-^o] + 6.3 1L24^ 1t1 [Lq]. 



(4) 



For the derivation of this relation, Zhu et al. (2008) relied on a 
sample of 187 star-forming (non-AGN) galaxies. 

Another relation, reported in Leroy et al. (2008), also com- 
putes the SFR from both FUV and 24yum data. Since a large frac- 
tion of our sample galaxies are characterized by a much higher star 
formation activity with respect to the sample in Leroy et al. (2008), 
we will apply the attenuation correction as reported in Zhu et al. 
(2008) and rely on the SFR relation in Kennicutt et al. (2009): 

1-28; 



SFR = 0.88 x 10~ is L F 



(5) 



where the units of Lp UV , c( , n are in er g S Hz -1 , to derive a star for- 
mation rate from the extinction-corrected FUV data. This relation 
was derived for a sample of normal galaxies in the present-day uni- 
verse and most star-forming galaxies out to redshifts z ~ 1 . With an 
upper limit of log Ltir[L q ] ~ 1 1.9, the relation is not immediately 
representive for the most dust-obscured LIRGs or ULIRGs found 
in the present-day universe (but this is not applicable for galaxies 
in our non AGN-dominated sample). 

Figure 4 shows the correlation between Z.^ C ii] and the 
SFRp UV+ 24^m derived in this way. The best fitting line through all 
data points is given by: 

log SFR FU v + 24„ m = 0.983 log L lCll] - 40.012 (6) 

where the units of SFR FU v + 24 f , m and L[ Cll ] are in M Q yr -1 and erg s , 




(erg .-)] 



Figure 4. SFRFt/v+24um (attenuation correction and SFR relation from Zhu 
et al. (2008) and Kennicutt et al. (2009), respectively) as a function of £[Cn]> 
central panel: SFRf-(/y + 24,, m /4;rZ) 2 as a function of [Cn] flux, bottom panel: 
residuals of the SFR/? (/v+24/im-£[Cuj plot. In this residual plot, Hn/starburst 
and LINER galaxies are indicated as blue filled squares and red diamonds, 
respectively. 



respectively. The spread around the mean trend is represented in the 
bottom panel of Figure 4 and can be quantified by a lcr dispersion 
of 0.27 dex. The central panel of Figure 4 again confirms that the 
observed trend is not due to a distance bias in our sample, since the 
lcr spread around the mean trend (0.26 dex) does not differ in this 
distance independent relation. 

Relying on the relation in Zhu et al. (2008), we can quan- 
tify the fraction of dust-enshrouded and unobscured star formation. 
Figure 5 plots the ratio between the observed FUV luminosity and 
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Table 3. Star formation rates for the galaxies in our sample, derived from 
different SFR relations based on MIPS 24/jm data, or a combination of 
GALEX FUV with MIPS 24pm data. 



Log[SFR ZBu (M 



Figure 5. The ratio Lfuv / 6.31 L24^m as a function of the total SFR, as 
derived from formula 4 and 5. This ratio of luminosities represents the frac- 
tion of unobscured star formation. The factor a in the denominator is the 
proportionality factor adopted from Zhu et al. (2008). For a = 6.31, they 
found that the sum LFUV.obs + 6-31 i24//m is a good approximation of the 
attenuation corrected FUV flux. 



the scaled MIPS 24pm luminosity, for which the scaling coefficient 
a = 6.31 was adopted from Zhu et al. (2008), as a function of the 
SFR. This Lpu V -to-6.31L24/ 1 m ratio represents the unobscured frac- 
tion of star formation. Figure 5 illustrates that the majority of our 
sample galaxies have a dominant contribution to the 24pm emission 
from dust grains heated by a young stellar population, rather than 
direct FUV emission from young stars. This makes the 24pm or a 
combination of GALEX FUV data and an appropriate attenuation 
correction a more reliable star formation indicator for our sample 
than solely relying on UV data. One exception is NGC4278, for 
which the unobscured fraction of star formation dominates dust- 
enshrouded 24pm emission of young stellar objects. Furthermore, 
Figure 5 shows a hint for a decreasing trend in the unobscured frac- 
tion of star formation towards higher SFR. This trend can be ac- 
counted for by an increasing opacity for galaxies which form stars 
at a higher rate. 



4 DISCUSSION 

4.1 [Cn] as a SFR indicator 

In the previous section, we found that the star formation rate cor- 
relates well with the [Cn] luminosity. Table 3 summarizes the star 
formation rates obtained from the two different star formation rela- 
tions for all galaxies in our reduced sample. Table 4 summarizes the 
coefficients a and b for the best fitting line y = ax + b for each star 
formation tracer, the uncertainty on the slope a and intercept b and 
the lcr dispersion of the individual galaxies around this mean trend. 
For both star formation relations, the uncertainties on the slope and 
intercept of the best fitting line are small and the spread around 
the mean trend is narrow. The dispersion around the correlation is 
smaller for the SFR derived from a combination of GALEX FUV 
and 24pm data (0.27 dex), than when estimating the SFR from the 
single 24pm luminosities (0.30 dex) (see also the bottom panels of 
Figures 3 and 4). Considering that a combination of FUV and 24pm 
data traces the complete star formation activity, as opposed to the 
single 24pm data, which only traces the obscured fraction of star 
formation, we will use the SFR tracing both the dust-enshrouded 
and unobscured activity to calibrate the SFR relation. Since the dis- 



Name 


Rieke2009 


Zhu20()8+Kennicutt2009 




(M Q yr 1 ) 


(M Q yr-') 


(1) 


(2) 


(3) 


Cartwheel 


3.25 ± 1.14 


7.97 ± 1.63 


NGC 0520 


5.04 ± 1.77 


7.02 ± 1.43 


NGC 0625 


0.04 ± 0.01 


0.06 ± 0.01 


NGC 0660 


1.23 ±0.43 


1.70 ± 0.35 


NGC 0695 


42.37 ± 14.88 


62.63 ± 12.78 


NGC 0986 


3.87 ± 1.36 


5.39 ± 1.10 


UGC 02238 


13.04 ±4.58 


18.19 ± 3.72 


NGC 1156 


0.06 ± 0.02 


0.18 ±0.04 


NGC 1266 


2.40 ± 0.84 


3.34 ± 0.68 


NGC 1317 


0.18 ±0.06 


0.31 ± 0.06 


NGC 2388 


17.31 ± 6.08 


23.35 ± 4.79 


NGC 4041 


1.82 ±0.64 


2.62 ± 0.53 


NGC 4 189 


0.25 ± 0.09 


0.42 ± 0.09 


NGC 4278 


0.014 ± 0.005 


0.03 ± 0.01 


NGC 4293 


0.46 ±0.16 


0.64 ±0.13 


NGC 4299 


0.18 ±0.06 


0.52 ±0.11 


NGC 4490 


0.26 ± 0.09 


0.42 ± 0.08 


NGC 4651 


0.47 ±0.17 


0.87 ±0.18 


NGC 5713 


6.04 ±2.12 


8.73 ± 1.78 


NGC 6240 


101.85 ± 35.78 


127.15 ± 25.94 


NGC 7217 


0.30 ±0.11 


0.49 ±0.10 


NGC 7552 


10.99 ± 3.86 


15.26 ±3.11 


NGC 7714 


8.85 ± 3.11 


13.34 ± 2.72 


NGC 7771 


14.06 ± 4.94 


19.47 ± 3.97 



persion is smaller for the SFR when estimated from both the FUV 
and 24pm luminosity (0.27) rather than from the 24pm data (0.30 
dex), we believe this as a confirmation of the link between the [Cn] 
emission and star formation activity in a galaxy, that has been re- 
ported in this work. 

From Eq. 6 we derive the SFR calibration: 

(Ac,,]) - 983 



SFR : 



(7) 



1.028 x 10 40 

where the SFR and L[Cn] are in units of M Q yr -1 and erg s" 1 , respec- 
tively. The SFR calibration factors are derived assuming a Kroupa 
(2001) IMF. This relation is valid for star-forming, late-type galax- 
ies with a star formation activity in the range 0.03 - 127 M Q yr~' 
and objects with a [Cn] luminosity between 39.3 < log L[ C u][ er g 
s -1 ] < 42.7 and a TIR luminosity between 42.1 < log Lri R [erg s -1 ] 
<, 45.4. 

Since the residual plot of this SFR relation (see Figure 4, bot- 
tom) indicates that the LINER galaxies in are sample contribute 
most to the observed scatter, we investigate whether a possible 
AGN contribution to the LINER galaxies introduce a bias to the 
derived SFR relation. If we calibrate the SFR relation for the 16 Hn 
and starburst galaxies in our sample, we should see a significant 
difference if the LINER sample indeed introduces a bias, possibly 
due to a small AGN contribution to the 24pm luminosities. From 
Figure 6 we can deduce that the best fitting line in our sample does 
not change substantially if we neglect the LINER galaxies for the 
SFR calibration. This already shows that the LINER sample does 
not introduce any bias due to a small AGN contribution. Moreover, 
half of the LINER galaxies are situated below the best fitting line. 

From this tight correlation between the star formation rate and 
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Table 4. Coefficients and uncertainties of the best fitting line for the relation between L[Cu] and two SFR relations. 

Variable Slope Intercept l<r dispersion (dex) 



SFR 24 ^ m ,ii,>fe2009 1.040±0.035 -42.535± 1.441 0.30 

SFRf (/v+2V „,z;, I ,2008 + ^™,„«,2009 0.983 ±0.021 -40.012±0.858 0.27 



Hll/starburst 
LINER 

Best fit (Hll/starburst) 
Best fit (HII/starburst + LINER) 




L°9[L[cnl (erg s" 1 )] 



Figure 6. SFRFi/y+24um (attenuation correction and SFR relation from 
Zhu et al. 2008 and Kennicutt et al. 2009, respectively) as a function of 
£[Cu]- In this pl°l. Hn/starburst and LINER galaxies are indicated as blue 
filled squares and red diamonds, respectively. The best fitting line for the 
Hn/starburst sample and the combined Hn/starburst+LINER sample are 
shown as blue plain and red dashed-dotted lines, respectively. 



a this paper 

Boselli et ol.(2002) 



, o'o 



Lo9[L[c] (erg s")] 

Figure 7. The star formation rate as estimated from relation 8 in this paper 
(blue filled squares) and from the SFR calibration in Boselli et al. (2002). 
Since Boselli et al. (2002) assumed a Salpeter (1955) IMF (a = 2.35) in 
the mass range between 0.1 and 100 M G , we have divided their calibration 
coefficient by 1.51, to convert it to the Kroupa (2001) IMF we applied. The 
best fitting lines for our calibration and the one in Boselli et al. (2002) are 
indicated as blue dotted and red dashed lines, respectively. 



the [Cn] luminosity, we conclude that the [Cn] luminosity is a good 
star formation rate indicator. The applicability of [Cn] as a star for- 
mation rate diagnostic had already been hinted by several authors 
(Stacey et al. 1991; Pierini et al. 1999; Leech et al. 1999; Boselli 
et al. 2002; Stacey et al. 2010), but we were able to quantify this 
correlation in an accurate way. Figure 7 compares our SFR relation 
to the calibration obtained in Boselli et al. (2002). Since Boselli 
et al. (2002) assumed a Salpeter (1955) IMF (a = -2.35) in the 
mass range between 0.1 and 100 M Q , we have divided their cali- 
bration coefficient by 1.51, to convert it to the Kroupa (2001) IMF 
(see Calzetti et al. (2010) for the derivation of the factor 1.51 dif- 
ference between the Salpeter (1955) and Kroupa (2001) IMF cali- 
brations). For low [Cn] luminosity objects both calibrations are still 
quite consistent, but for galaxies with an increasing [Cn] luminosity 
the SFR estimate provided by Boselli et al. (2002) quickly diverges 
from our estimates up to a factor of ~5 for the highest luminosity 
objects in our sample. This deviation towards higher [Cn] lumi- 
nosities is probably due to an underestimation of the Her extinction 
in more luminous galaxies, since the attenuation relations at that 
time where calibrated without taking into account heavily obscured 
star formation and thus underestimating the true Ho- emission. The 
scatter around the mean trend in our SFR relation is smaller than 
a factor of 2, while the dispersion in the [Cn]-H„ luminosity re- 
lation already reaches a factor of ~4. The final uncertainty on the 
SFR calibration in Boselli et al. (2002) is estimated to be as high 
as a factor of ~ 10. We believe the poorly known characteristics of 
the [Nn] contamination and the attenuation correction at the time 



Boselli et al. (2002) performed their SFR calibration explain a sig- 
nificant amount of their reported scatter. Benefiting from the in- 
creased database of [Cn] observations and exploring other reliable 
SFR tracers in this paper, we can revise the analysis in Boselli et 
al. (2002) and conclude that the [Cn] emission is a reliable SFR 
indicator in most normal star-forming galaxies. 

4.2 Applicability 

The extension of this SFR relation to ULIRGS should be treated 
with caution, since a decrease of the Lr Cll ]-to-L F T R ratio with increas- 
ing warm infrared color for all galaxy types has been observed in 
many samples (Crawford et al. 1985; Stacey et al. 1991; Malhotra 
et al. 1997, 2001; Luhman et al. 2003; Verma et al. 2005; Brauher et 
al. 2008). More specifically, Luhman et al. (2003) report a [Cn] line 
deficit in a sample of 15 ULIRGS, after they had noticed this trend 
in a few individual ultra-luminous objects (Luhman et al. 1998). 
Figure 8 hints at a similar, but weaker, trend for the galaxies in our 
sample. This effect is probably due to a decrease in efficiency of the 
photoelectric heating of the gas in strong radiation fields (Malhotra 
et al. 2001). Other possible explanations for this trend were given 
in Negishi et al. (2001) and Bergvall et al. (2000). Negishi et al. 
(2001) suggest an increased collisional de-excitation of [Cn] due to 
an enhanced gas density, or a decrease in the ionized component 
for increasing star formation activity. While Bergvall et al. (2000) 
invokes the self-absorption of [Cn] for galaxies with increasing 
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Figure 8. The ratio LfCuj-to-LpiR as a function of the effective dust tem- 
perature, as embodied by the / r (60yum)-to-/,,(100 J um) ratio, where the FIR 
luminosity was derived following Helou et al. (1985). We notice a decrease 
of the L[c u ]-to- Lfir ratio with increasing warm dust temperature. 

metallicity. Recently, Papadopoulos et al. (2010) and Rangwala et 
al. (2010) claimed that a high dust optical depth could be the cause 
of this [Cn] line deficit in several galaxies. If this trend in the [Cn] 
line deficit holds for more IR-luminous galaxies, this might have its 
implications for the reliability of [Cn] as a star formation indicator 
in those objects. Even more recently, Gracia-Carpio et al. (2011) 
inferred from their Herschel observations a similar line deficit for 
several other FIR fine structure lines ([Nn], [Oi] , [Nm] , [Om]). 
They allocate this deficit to a transition between two modes of star 
formation (in normal disk galaxies and major merger systems) with 
a different star formation efficiency. 

On the other hand, Pierini et al. (1999) also report a decrease 
in the L[ Cl ,]-to-L F i R ratio for decreasing Ho- EW. This trend has been 
observed for more quiescent galaxies with Ho- EW < 10 A and 
can be explained by an increasing contribution from older, less 
massive stars that heat dust grains and therefore contribute to the 
FIR-luminosity. The radiation from these stars is not hard enough 
to substantially heat the dust grains through the photoelectric ef- 
fect (Pierini et al. 2003). Since collisions with these photoelectrons 
are the main heating source for the interstellar gas, this inefficient 
gas heating will render those galaxies [Cn] quiet for an increas- 
ing contribution from a more evolved stellar population. This effect 
makes the extension from the SFR relation to quiescent galaxies 
not straightforward either. For objects with an extremely low star 
formation activity, Pierini et al. (1999) also found that the main 
contribution to the [Cn] emission arises from the diffuse neutral in- 
terstellar medium. In our current sample, two galaxies (NGC 4293 
and NGC 7217) satisfy the criterion Ha EW < 10 A, characteriz- 
ing a quiescent galaxy (Kennicutt & Kent 1983; Kennicutt et al. 
2008). Since those galaxies do not deviate from the mean trend in 
the SFR-L[ Cll j plot, a larger sample of quiescent galaxies is neces- 
sary to further investigate whether our SFR relation is expandable 
to objects with a lower SFR compared to the galaxies in the present 
sample. 

Because of systematic variations in the Lfdj-to-LpiR ratio 
and varying contribution with luminosity from several components 
(PDRs, diffuse cold neutral and warm ionized gas, Hn regions) to 
the [Cn] emission, we should be careful in extending the SFR re- 
lation in this paper to more extreme luminosities, on both ends of 
the luminosity scale. Although the derived star formation relation 
might not be valid for more extreme cases of star formation activ- 



ity (either very quiescent or actively star-forming), this calibration 
might be useful as an alternative indicator of the star formation 
rate in objects which lack either UV or IR observations. Since [Cn] 
emission is hardly not affected by attenuation in most cases (cur- 
rently, Arp 220 is the only object where prominent obscuration ef- 
fects have been claimed), it might provide an immediate probe of 
the complete star formation activity. 

At high redshift, the behaviour of the L[c u ]-to-LpiR ratio is even 
different. Ultra luminous objects (ULIRGs; Ltir > 10 12 L Q ) do not 
seem te be affected by a [Cn] deficit. On the contrary, Maiolino 
et al. (2009) suggest an enhancement by at least one order of 
magnitude of the [Cn] emission in galaxies at high-redshift com- 
pared to local galaxies of the same infrared luminosity. Maiolino 
et al. (2009) argues that this enhanced [Cn] emission could be due 
to lower metallicities of the ISM in those high-redshift galaxies, 
which tend to have a lower dust content and therefore a larger [Cn] 
emitting region. Independent of the physical origin of this effect, 
if confirmed by future observations, the strong [Cn] emission from 
IR luminous high-redshift objects would have great implications 
for future observations at high redshift. This suggested increase in 
detectability of the [Cn] line at high redshift would also imply an 
extensive applicability of a SFR relation based on [Cn] luminosi- 
ties, which we will derive in this paper. 

Recently, Stacey et al. (2010) reported the detection of the 
[Cn] line in 12 galaxies at redshifts ranging from 1 to 2. They 
concluded that for starburst-powered galaxies the [Cn] emission is 
comparable to values found in local star-forming galaxies with sim- 
ilar FIR luminosities ([Cn]/FIR ~ 3xl0~ 3 ), while AGN-dominated 
galaxies are characterized by [Cn]-to-FIR luminosity ratios similar 
to local ULIRGs ([Cn]/FIR ~ 4x 10"), suggesting that this effect of 
enhanced [Cn] emission only becomes apparent at even higher red- 
shift (z>2). Stacey et al. (2010) also reported that [Cn] is a reliable 
star formation indicator for their heterogeneous sample of starburst- 
and/or AGN-dominated galaxies in the redshift interval z = [1,2]. 
At even higher redshifts (z>2), the lower metallicity might influ- 
ence the temperature and chemical structure of PDRs and thus, the 
[Cn] emission (e.g. Lequeux et al. 1994; Wolfire 1995; Rollig et al. 
2006). Future observations of ULIRGs with Herschel and ALMA 
will provide insight in the applicability of [Cn] as a star formation 
indicator at high redshift and the validity of the star formation rela- 
tion derived in this paper. 

4.3 Nature of the [Cn] emission in galaxies 

Besides the first accurate quantitative calibration of the star forma- 
tion rate against the [Cn] luminosity, the tightness of this correla- 
tion gives us insight in the origin of the [Cn] emission on a global 
galaxy-scale. In this section, we introduce two possible explana- 
tions for the tightness of the SFR-[Cn] relation. 

Although PDRs are the main contributor to the [Cn] emission 
in most galaxies, a significant fraction also originates in the cold 
neutral medium (CNM) (i.e. Hi clouds, Kulkarni & Heiles (1987); 
Madden et al. (1993); Bennett et al. (1994); Wolfire (1995); Sauty 
et al. (1998); Pierini et al. (1999, 2001); Contursi et al. (2002)), 
the warm ionized medium (WIM) (i.e. diffuse Hn regions, Madden 
et al. (1993); Heiles (1994); Malhotra et al. (1997, 2001); Leech 
et al. (1999); Contursi et al. (2002)) and to lesser extend also in 
Hn regions (Stasinska 1990) (see Boselli et al. 2002 for a quantita- 
tive analysis of these different contributions). Since the contribution 
from compact Hn regions to the [Cn] emission is negligible with re- 
spect to other components on galactic scales (e.g. Stasinska 1990; 
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Malhotra et al. 2000), it is difficult to believe that the contribution 
from these star-forming regions alone causes this good correlation. 

This brings us to a first possible explanation for the strong 
correlation between [Cn] and the SFR. Considering that PDRs are 
neutral regions of warm dense gas at the boundaries between Hn 
regions and molecular clouds, we think that most of the [Cn] emis- 
sion from PDRs arises from the immediate surroundings of star- 
forming regions. This might not be surprising, since FUV photons 
from young O and B stars escaping from the dense Hn regions, 
impinge on the surface of these neutral PDR regions where they 
will heat the gas through the photoelectric effect on dust grains. 
We believe that a more or less constant contribution from PDRs to 
the [Cn] emission and the fact that this [Cn] emission from PDRs 
stems from the outer layers of photon-dominated molecular clumps 
(i.e. at the boundary of molecular clouds and Hn regions) might be 
a reasonable explanation for the tight correlation between the star 
formation rate and L[Cn]- Also, Malhotra et al. (2001) suggest this 
tight correlation between PDRs and star-forming regions from their 
analysis of the FUV flux Go and gas density n in PDRs. The high 
PDR temperature and pressure required to fit their data, imply that 
most of the line and continuum FIR emission arises from the im- 
mediate proximity of expanding Hn regions. Moreover, mapping 
of [Cn] in the Milky Way (Bennett et al. 1994) and in spatially re- 
solved nearby galaxies (see Sauvage et al. 2005 for an overview) 
gives indications for a close association between PDRs and ionized 
gas in these galaxies. 

Alternatively, we consider the likely possibility that this tight 
correlation is not the reflection of shared photoexcitation processes, 
taking place at the same position within a galaxy. On the contrary, 
the [Cn] emission might not be directly linked to the star forma- 
tion, but could instead trace the cold ISM being therefore indi- 
rectly related to the star formation through the Schmidt law (Kenni- 
cutt 1998). Evidence confirming this indirect link is the association 
of [Cn] emission with modest densities and softer FUV radiation 
fields (e.g Pineda et al. 2010). 

A similar connection between the [Cn] emission in a galaxy 
and the diffuse ISM is typically present in more quiescent objects 
(Pierini et al. 1999, 2001, 2003), while in star forming galaxies the 
[Cn] line emission is found to mainly arise from PDR's (Crawford 
et al. 1985; Stacey et al. 1991; Madden et al. 1993). Mapping of 
[Cn] at high resolution is necessary to ascertain the nature of the 
tight correlation between the SFR and [Cn] emission for the normal 
star-forming galaxies in our sample. 

Future observations of spatially resolved objects in the nearby 
universe and at high redshift of [Cn], CO(1-0) and [Nn] lines for 
a sample of galaxies ranging several orders of magnitude in IR- 
luminosity, will be able to disentangle the different components 
contributing to the [Cn] emission, since the CO(1-0) intensity cor- 
relates well with [Cn] (Hollenbach & Tielens 1999) in PDRs and 
[Nn] traces the Hn regions and the warm ionized medium (Wright 
et al. 1991; Bennett et al. 1994). Such analysis will enable us to 
make a distinction between the different sources that contribute to 
the [Cn] emission and, in particular, examine the tight correlation 
between the star formation rate and [Cn] luminosity. High spatial 
resolution data will be able to distinguish whether [Cn] emission 
from PDRs is closely related to star-forming regions in galaxies 
or rather the global [Cn] emission in a galaxy traces the gas mass. 
The latter would imply an indirect link between [Cn] and the star 
formation activity in a galaxy through the Schmidt law. 



5 CONCLUSIONS 

The aim of this paper was to investigate the reliability of [Cn] as a 
SFR indicator. The analysis was conducted in the following way: 

• We assembled a sample of 38 galaxies for which ISO [Cn], 
GALEX FUV and MIPS 24 /jm data are available. This initial sam- 
ple was subdivided according to their nuclear spectral classifica- 
tion (Hn/starburst, LINER, AGN). Due to the contaminating 24yum 
emission in AGN hosts, we were unable to directly relate the 24yum 
emission in those objects to the star formation activity. Therefore, 
the objects with nuclear spectra resembling those of an AGN, were 
eliminated from our sample. 

• For the remaining 24 galaxies in our sample, the SFR was es- 
timated from either the 24/im luminosity or a combination of FUV 
and 24yum data. In the latter case the FUV emission is corrected for 
internal dust attenuation according to Zhu et al. (2008). The SFR 
calibration in Kennicutt et al. (2009) provides a SFR estimate for 
this extinction-corrected FUV luminosity. Comparing the disper- 
sions for these two different SFR relations, we found that the SFR 
correlates best with the [Cn] luminosity, when the star formation 
activity is traced by a combination of those luminosities. 

• From the tight correlation of the SFR with the [Cn] luminosity 
(the dispersion around the mean trend is 0.27 dex), we conclude 
that the [Cn] luminosity is a reliable SFR diagnostic in normal, 
star-forming galaxies in the local universe and find the following 
SFR calibration: 

(Wergs-']) 0983 
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• The extension of this relation to more quiescent galaxies (Ha 
EW < 10 A) and ultra luminous galaxies (ULIRGs) should be 
treated with caution, due to a deviation from the L[ Cll ]-to-L F i R cor- 
relation at those extreme ends of the luminosity scale. 

• For galaxies in the redshift range z= [ 1 ,2] , the comparable [Cn] 
luminosities and the enhanced [Cn] emission, at even higher red- 
shift, compared to local galaxies with similar FIR luminosities, will 
render an immense applicability of this SFR relation in future high 
redshift surveys (both with Herschel and ALMA). 

• We believe this tight correlation between the [Cn] luminosity 
and the star formation activity may imply a direct link in photo- 
excitation processes: the [Cn] emission from PDRs arises from the 
immediate surroundings of star-forming regions (i.e. from the outer 
layers of photon-dominated clumps) and contributes a more or less 
constant fraction on a global galaxy-scale. Alternatively, the [Cn] 
emission might equally well trace the star formation in a galaxy, 
when assuming that the [Cn] emission is associated to the cold ISM 
in a galaxy and, therefore, is indirectly linked to the SFR through 
the Schmidt law Esfr = 
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